Crystal Growth as an Excitable Medium 
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Abstract 

Crystal growth has been widely studied for many years, and, since the pioneer- 
ing work of Burton, Cabrera, and Frank, spirals and target patterns on the crystal 
surface have been understood as forms of tangential crystal growth mediated by 
defects and by two-dimensional nucleation. Similar spirals and target patterns are 
ubiquitous in physical systems describable as excitable media. Here we demon- 
strate that this is not merely a superficial resemblance; that the physics of crystal 
growth can be set within the framework of an excitable medium, and that appreciat- 
ing this correspondence may prove useful to both fields. Apart from solid crystals, 
we discuss how our model applies to the biomaterial nacre, formed by layer growth 
of a biological liquid crystal. 

1 Introduction 

Molecular- scale spiral and target patterns on the surface of crystals were predicted 
even before they were seen directly with microscopes. Burton, Cabrera, and Frank 
proposed a dynamical model to explain how crystals should grow layer by layer 
through the addition of material at kinks; growth sites at the edges of a spiral or 
target pattern 1 1 1. The Burton-Cabrera-Frank (BCF) model has since been am- 
ply confirmed by observations, although today the quantitative dynamics of crystal 
growth is known to be affected by a host of effects not considered in the minimal 
BCF model. On the other hand, spiral and target patterns similar to those on the 
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surface of crystals are seen in chemistry, e.g., in the Belousov-Zhabotinsky reac- 
tion; in biology, e.g., in yeast growth or in heart fibrillation; and in many other in- 
stances. The latter phenomena are all examples of excitable media, a rather general 
class of systems in which elements are quiescent until excited by some stimulus, 
after which they are unresponsive to further stimuli during some refractory period 
before returning to their initial quiescent, excitable state (2][3]. This very sim- 
ple physics is sufficient to produce complex spatiotemporal patterns of targets and 
spirals, which is why these patterns are observed so generally in many different 
fields [4,5,6,7,8]. The resemblance of spiral and target patterns in crystal growth 
to such patterns is not mere accident. Here we show with a coupled map lattice 
model that the classical BCF model of crystal growth is also an excitable medium. 
As well as tangential growth of a solid crystal, we demonstrate that the excitable 
paradigm holds and the BCF model may be applied to layer by layer growth of a 
liquid crystal, which has recently been noted in the biomaterial nacre (9). 

2 The BCF model 

Crystals grow either layer by layer, in tangential growth, or continuously, in normal 
growth (10]|TT). In tangential growth the crystal surface is smooth on a molecular 
scale and consists of crystallographic planes of layers of growth units (atoms, ions, 
or molecules depending on the crystalline material). Burton, Cabrera, and Frank 
described how incomplete layers are bounded by steps, and incomplete rows of 
growth units along a step end at kinks. New growth units are incorporated at kinks. 
New steps are created only by two-dimensional nucleation: the creation of a new 
island layer on top of an existing layer. The other mode of growth of such a layer 
structure is by a screw dislocation on the crystal face. In this case, a step propagates 
upwards in a helicoidal pattern, like a spiral staircase. A general point on the 
surface is then part of an already formed lower layer, and is quiescent, but receptive 
to the nucleation of a new layer above it given a large enough perturbation, that is 
to say a fluctuation that concentrates sufficient material nearby to nucleate a new 
surface. The edges of new layers are where new material is generally incorporated, 
while just after a new layer has formed at some point there is a refractory period 
during which another new layer cannot be formed, because any new material will 
of preference be incorporated at the growing edge nearby. 

3 Coupled map lattice model 

Layer growth proceeds by the incorporation of individual growth units — it is es- 
sentially a discrete process — so we employ a discrete model. (An early argument 
for introducing such discrete models into crystallography was given by Mackay 
1 13].) Thus we construct a minimal coupled map lattice model of the growth dy- 
namics (Fig. [TJ. The surface is divided up into cells, which may be thought of as 
the size of the growth units that compose the crystal. In order to avoid anisotropic 
growth, we use a randomized grid, as introduced by Markus and Hess 1 14], and, 
defining a neighbourhood radius, R, the neighbours are all elements within this 
radius 

\rij -r k i\ < R, 
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where ij and kl are the cell coordinates on the bidimensional lattice and we use the 
length of one cell as the length unit. The parameter R will influence the separation 
between terraces on the growing crystal surface, and so physically we can associate 
it with the surface diffusion of the depositing particles. Each cell has an associated 
height, Hij , initially zero and updated at the end of each time iteration. Hij is 
a continuous variable into which a random component is introduced, as we show 
below, since there has to be certain variability for defects to occur. Cells can be 
considered to be in an excitable state or in a refractory state. The excitable state 
means that the cell may nucleate a new island or add new material at a growth front. 
The condition for nucleation is that the cell in question must be on a flat surface, 
meaning that the height difference with its neighbouring cells must be smaller than 
a certain margin, AHn'. 

AH<AH N . 

neighbors 

This parameter will define the frequency of nucleations and target patterns. In 
crystallization it is related with the supersaturation or supercooling of the precipi- 
tating material. The condition for growth is that the cell must be at the edge of a 
growth front, meaning it has at least one neighbour with a height difference larger 
than a certain threshold, AHq : 

AH > AH G . 

This parameter is responsible for the appearance of screw dislocations. Physically 
there can be several causes for this kind of defects, such as the flexibility of the 
bonds of the lattice, the heterogeneity of growth units (for example in protein- and 
virus-crystal growth, and in liquid crystallization), and the presence of impurities. 
If the cell does not fulfil either of these conditions, then it is considered to be in 
a refractory state. The duration of this state can be altered by increasing or de- 
creasing the radius of the neighbourhood; a larger radius implies more neighbours 
and hence a longer refractory period, which is reflected in the separation between 
terraces in the spiral and target patterns. However, a larger neighbour radius would 
also imply a higher probability of being at the edge of a step. This is reflected 
in the front thickness, meaning that the number of cells that experience growth in 
each iteration is greater, giving the appearance that the front is spreading faster. 
But since our time units (iterations) are completely arbitrary, the kinetics of crystal 
growth are beyond the focus of our model. 

If the cell is in a position to nucleate, it must first pass a probability check 
with probability Pat, as nucleation is a stochastic process. Pn is generally small 
(<C 1). If the check is positive and nucleation takes place, the height of the cell is 
increased by 1 plus or minus a small random factor a: 

H ij (t + l) = H ij (t) + l±a. 

On the other hand, when growth takes place, the height of the cell is increased by 
the mean height difference with its higher neighbours: 

Hij (t + 1) = H l3 (t) + £ H ^)- H ^) ± a , 

k,i 71 

where (A;, I) are the coordinates for each higher neighbour and n is the total number 
of higher neighbours, which depends on R. This growth algorithm is performed 
simultaneously for all cells, and the process is iterated. 
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Table 1 : Model parameters 





Description 


Influence 


Physical meaning 


a 


Stochastic term 


Lattice defects 


Impurities 


AH N 


Height margin to consider flat 


Nucleations and target pat- 


Supersaturation 




enough to nucleate 


terns 




AH G 


Height threshold to consider the 


Screw dislocations and spi- 


Particle shape variability 




base of a step 


ral patterns 




R 


Radius of the cell neighbourhood 


Terrace separation 


Surface diffusion 



4 Results 

When a nucleation event occurs it automatically inhibits further nucleation in its 
neighbourhood and all new material will be added at the edge of the newly formed 
island. If this island can grow through several time steps without greatly varying its 
height, a second island may nucleate on top of it. If this process continues period- 
ically it will produce a target pattern (Fig. [2^). In order to produce this pattern we 
set a broad nucleation margin A if at to favour nucleation and a low stochastic term 
a to avoid rough surfaces. In crystallization these parameters would respectively 
translate to a high supersaturation or supercooling of the precipitating material 
and a low shape variability of its growth units. Although we can currently only 
describe the qualitative relation between the model and the physical parameters, 
understanding this behaviour is necessary in order to obtain a quantitative relation. 

Another common pattern in both excitable media and crystal growth is the 
spiral (Fig. [2|d). When two growth fronts with slightly different heights collide, 
a portion of one may overlap the other, and continue its growth revolving around 
the dislocation centre. These patterns are produced with a small active neighbour 
threshold AHg, which allows growth fronts to overlap instead of annihilating each 
other, and a larger stochastic term a that introduces variability in the front height. 
Physically these parameters would imply a crystalline lattice with high tendency 
to create dislocations of whatever type, whether through the presence of impurities 
or due to the size and shape variability of the growth units. 

An important aspect of the behaviour of excitable media is that excitable waves 
cannot pass through each other, but rather destroy each other, as the refractory 
period does not permit wave propagation to continue; nor do waves reflect from 
boundaries. In the well-known excitable model of forest fires, for example, a sec- 
ond fire cannot pass until the vegetation burnt by the first has grown back, and 
excitable cardiac cells cannot fire again until they have recovered from their earlier 
firing. In Fig. [3^ we see a sequence of two growth fronts colliding and annihilating 
each other. In this case the height of both fronts is very similar, as the stochas- 
tic term a has been set to be small, and the threshold for overlapping, AHg, is 
too large for one front to grow over the other. With a higher variability in the 
front height and imposing a less restrictive growth threshold, one of the fronts 
may absorb the other, producing an edge dislocation, or it may partially overlap it, 
producing a screw dislocation (Fig. [3]}). 

An unusual aspect of our growth model compared to other excitable media is 
that here the third dimension of space — outwards from the material — is equiv- 
alent to time. Thus one can see in three dimensions the trajectory of the core of 
a spiral — a screw dislocation in the material. We can follow the position of this 
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core by calculating the Burgers vector [15], the anholonomy in the lattice. The 
traditional manner of calculating this vector is by drawing a Burgers circuit around 
the defect and comparing it to the same circuit in a perfect lattice. But we cannot 
do this in our case because the lattice is not perfectly periodic: it has a random 
component on the plane and a continuous coordinate, H, in the third dimension. 
Instead we take the following approach: We follow a square circuit around each 
cell looking for large height discontinuities (> AHg) that indicate the presence 
of a step. We accumulate the height differences around the whole circuit, adding 
when going up a step and subtracting when going down, and the resulting number 
is the magnitude of the Burgers vector in the vertical direction. If there is no defect 
around the considered cell, then the Burgers vector should be zero. But, if we are at 
the centre of a screw dislocation, the magnitude of the vector will be different from 
zero and its sign will mark the sense in which the spiral rotates (Figs [4^ and b). 
As the spiral continues growing, the position of its centre oscillates back and forth 
following a roughly straight trajectory (Fig. [4J:). When we plot the distance from 
the spiral core to the origin of our map in this case, we can see that the oscillation 
is periodic and we can estimate the amplitude (~12 cells) and frequency (~1 cycle 
per 20 time units) of the oscillation (Fig. |4ji). In Figs[4|e, f) we can see another 
example where we use the Burgers circuit to find the two cores of a double spiral; 
a so-called Frank-Read source 1 16 1. The sign of the Burgers vector at the spiral 
cores is related to the rotation sense; a negative vector means clockwise rotation, 
and positive anticlockwise. 

5 Discussion 

Tangential crystal growth involving spirals and target patterns has been intensively 
studied in the decades since Burton, Cabrera and Frank's pioneering work 1 17 1. 
Falo et al. 1 18 1 used a Langevin approach; a sine-Gordon equation plus noise to 
simulate spiral growth plus nucleation. Aranson et al. |19 | pursued a Ginzburg- 
Landau approach for spiral growth. And Smereka [20] took a spiral-growth con- 
tinuum approach to BCF. More recent efforts include the application of phase field 
models that allow the quantitative calculation of physical properties such as kinetic 
coefficients [21 1 and the computation of more complicated morphologies such as 
dendrites 1221 . 

In contrast with the above-cited works, our discrete model puts both growth 
modes, of spiral growth and island nucleation, on an equal footing. In other words, 
we use a coupled map lattice model instead of partial differential equations because 
we can produce both spiral and target patterns without requiring any other symme- 
try breaking elements. The connection that we demonstrate here with the theory 
of excitable media has not previously been made. The connection is potentially of 
great utility to both fields, as excitable media have a series of theoretical models 
that can provide useful insight into crystal growth processes, while crystal growth 
is studied in much quantitative detail and can provide data with which to study it as 
an excitable medium. It is true that there are important crystal-growth phenomena 
that are not present in the minimal excitable-medium model; the Schwoebel effect 
leading to step bunching, etc [23, 24]. These act to make some crystal surfaces 
behave rather differently from the simple model. However many crystal surfaces 
do show the surface morphologies that reveal that this excitable dynamics is act- 
ing, as we may appreciate from scanning electron microscope and atomic force 
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microscope images, such as those of macromolecular crystals ifTOl ITTI WTil ; see 

The possibility that liquid crystals might grow layer by layer in the manner 
of a solid crystal seems to have been unexploited by man; not however by na- 
ture. In fact our interest in an excitable-medium description of layer growth arose 
out of our work with the biomaterial nacre l25l . Nacre, or mother of pearl, is 
the iridescent material that makes up the inner layer of many mollusc shells, and 
also the pearls that form around a foreign body within the mollusc. It is a mate- 
rial that has been much studied, because of its beauty, its strength, and its interest 
for biomimetics as a complex material self assembled outside the cells of the or- 
ganism that produces it. We have recently shown that nacre construction begins 
with the self assembly into a cholesteric liquid crystal of rod-shaped crystallites of 
the polysaccharide /3-chitin (9). Liquid crystals are usually formed in laboratory 
experiments by altering a global experimental variable, be it temperature, concen- 
tration, electric field, etc, that causes the whole of the experimental domain to 
crystallize at once. But the mollusc uses another mode of construction, in which 
the liquid crystal is built up layer by layer, in the manner of a solid crystal. It is 
for this reason that we find in nacre (see Fig. [6) the characteristic spiral, target, and 
labyrinthine patterns we have discussed above in the context of solid crystals. It 
is our opinion that understanding this equivalence between layer growth in solid 
and liquid crystals — and it appears probable to us that nature uses layer growth 
of liquid crystals ubiquitously — might provide useful new ideas for technological 
applications of liquid crystals. 

We would like to thank Carlos Pina for his permission to use the AFM im- 
ages of Figure 5. We acknowledge MINCINN (Spain) projects FIS20 10-22322- 
528C02-02 and CGL2010-20748-C02-01. 
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X AH<AH N 
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nucleation 



H ij (t+l) = H ij (t) + l±a 



growth 



AH > AH~ 




Figure 1: Flowchart of our coupled map lattice model combining the principles of 
tangential crystal growth and excitable media. 




(a) 



(b) 



Figure 2: (a) Target pattern (a = 0.08, AH N = 0.1, P N = 10~ 3 , AH G = 0.8, 
R = 10). (b) Spiral (a = 0.05, AH N = 0.06, P N = 1 x 10" 6 , AH G = 0.8, R = 6). 
The colour scheme represents height. 
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Figure 3: (a) Two growth fronts collide and annihilate each other (a = 0.05, AHjy = 
0.05, P N = 10~ 6 , AH G = 0.8, R = 6). (b) Partial overlap of two growth fronts 
produces a screw dislocation (a = 0.05, AH N = 0.05, P N = 10" 6 , AH G = 0.6, 
R = 6). 




Figure 4: (a) Spiral created by a screw dislocation (a = 0.05, AH^ = 0.05, Pn = 
10 -7 , AHq = 0.6, R = 5). (b) Burgers vector calculated on the same map. The value 
of the vector is null everywhere except at the centre of the spiral, (c) Trajectory of the 
centre of the spiral during 150 iterations, (d) Distance from the origin to the centre of 
the spiral versus time during 150 iterations, (e) Double spiral pattern created by screw 
dislocations (a = 0.05, AH N = 0.05, P N = 10~ 6 , AH G = 0.6, R = 4). (f) Burgers 
vector calculated on the same map. The sign of the vector at the spiral cores is defined 
by the rotation sense. 
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Figure 5: Spirals and target patterns in crystals: (a) spiral pattern on the surface of 
anglesite (PbS04) growing on a celestine (SrS04) crystal; scan area is 10 x 10 /mi ; 
(b) target pattern on the surface of a calcite (CaCOa) crystal growing in the presence of 
magnesium; scan area is 6 x 6 /im. Both images were taken during growth with AFM 
in contact mode lfT2ll . 



Figure 6: Spirals and target patterns in nacre; scanning electron micrographs of the 
growth surface of bivalve nacre: (a) Spiral (I so gnomon legumen) (b) Paired spirals 
{Pteria avicula); (c) Target pattern (Pteria avicula). 
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